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Abstract. Lithium abundances and rotation, determined 
for 120 subgiant stars in Lebre et al. (1999) are analyzed. 
To this purpose, the evolutionary status of the sample as 
well as the individual masses have been determined using the 
HIPPARCOS trigonometric parallax measurements to locate 
very precisely our sample stars in the HR diagram. We look 
at the distributions of A^and Vsini with mass when stars 
evolve from the main sequence to the subgiant branch. 

For most of the stars in our sample we find good agree- 
ment with the dilution predictions. However, the more mas- 
sive cool stars with upper limits of Li abundances show a 
significant discrepancy with the theoretical predictions, even 
if the Non-LTE effects are taken into account. For the ro- 
tation behaviour, our analysis confirms that low mass stars 
leave the main sequence with a low rotational rate, while 
more massive stars are slowed down only when reaching the 
subgiant branch. We also checked the connection between 
the observed rotation behaviour and the magnetic braking 
due to the deepening of the convective envelope. Our results 
shed new light on the lithium and rotation discontinuities in 
the evolved phase. 

Key words: convection - stars: abundances - stars: evo- 
lution - stars: interiors - stars: late-type - stars: rotation 



1. Introduction 

This paper is the second in a series about the study 
of lithium and rotation in evolved stars based on both 
new high resolution spectroscopic observations and pre- 
cise rotational velocities obtained with the CORAVEL 
spectrometer. In Lebre et al. (1999, hereafter Paper I) we 
derived lithium abundances by spectral synthesis analy- 
sis and presented the observational data for a large and 
homogeneous sample of F, G, and K-type Population I 
subgiant stars. On the basis of this uniform data set, we 
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could confirm the occurrence of the rotational disconti- 
nuity near the spectral type F8IV (Gray & Nagar 1985, 
De Medeiros & Mayor 1989, 1990), and localize a lithium 
drop-off around the spectral type G2IV. No clear corre- 
lation appeared between the lithium abundance and the 
rotational velocity (see also De Medeiros et al. 1997). 

In the present work we investigate the physical pro- 
cesses that underline the lithium and rotational disconti- 
nuities along the subgiant branch. For this purpose we first 
determine the evolutionary status and individual masses 
of our sample stars by using the HIPPARCOS parallaxes 
and by comparing the observational Hertzsprung-Russell 
diagram with evolutionary tracks computed with the Tou- 
louse-Geneva code (§2). This allows us to study very pre- 
cisely the behaviour of the lithium abundance and of the 
rotational velocity as a function of effective temperature, 
stellar mass, metallicity and evolutionary stage. The li- 
thium main features are discussed in §3 where we compare 
the observations with theoretical predictions related to 
the dilution mechanism (Iben 1967a, b). Finally, the con- 
nection between the observed rotation behaviour and the 
magnetic braking due to the deepening of the convective 
envelope (Gray 1981, Gray & Nagar 1985) is quantified in 
§4. This study, based on a close examination of the stellar 
parameters, sheds new light on the question of the link 
between rotation and lithium discontinuities in subgiant 
stars. 

2. Observational data and evolutionary status 

2.1. Working sample 

Our analysis is based on the observational data from Pa- 
per I. This sample is composed of 120 Pop I subgiant stars 
with F, G and K spectral types which belong to the "Cat- 
alogue of rotational and radial velocities for evolved stars" 
(De Medeiros & Mayor 1999), as well as to the Bright Star 
Catalogue (Hoffleit and Jaschek 1982). We thus use the ro- 
tational velocities given in Paper I, as well as the values 
derived for log g , and T c g with their respective errors. 
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Fig. 1. Distribution of the sample stars in the HR dia- 
gram. Luminosities and related errors have been derived 
from the Hipparcos parallaxes. The typical error on T e g is 
± 200 K (Paper I). Evolutionary tracks at [Fe/H]=0 are 
shown for stellar masses between 1 and 4 M Q . The turnoff 
and the beginning of the ascent on the red giant branch 
are indicated by the dashed and dotted lines respectively 
in order to discriminate dwarfs, subgiants and giants 

2.2. Evolutionary status 

In order to interpret accurately the observations, we need 
to know the mass and the evolutionary stage of the sample 
stars. We use the HIPPARCOS (ESA 1997) trigonometric 
parallax measurements to locate precisely our objects in 
the HR diagram. Among our 120 stars, only one object 
(HD 144071) has no available Hipparcos parallax and is 
thus rejected from further analysis. Intrinsic absolute mag- 
nitudes My are derived from the parallaxes and the my 
magnitudes given by Hipparcos. We determine the bolo- 
metric corrections BC by using the Buser and Kurucz's 
relation (1992) between BC and V-I (again taken from 
the Hipparcos Catalogue). Finally, we compute the stellar 
luminosity and the associated error from the sigma er- 
ror on the parallax. The uncertainties in luminosity lower 
than ±0.1 have an influence of ±0.4 in the determina- 
tion of the masses. We show the results of these determi- 
nations in Fig. and in Tab. [l| This table displays the 
stellar masses and luminosities derived for all objects of 
our sample. Moreover complementary data used for the 
present analysis (Vsini, T e g , log g , [Fe/H] and ) can 
be found in Papier I (see Tables 1 to 3). The error adopted 
on T c ff (±200K) is typical for this class of stars, as already 
discussed in Paper I. 

We have computed evolutionary tracks with the 
Toulouse- Geneva code for a range of stellar masses be- 




Log (T.„) 



Fig. 2. Distribution of the Li abundances in the HRD. 
Single and binary stars are identified with open and filled 
circles respectively. The size of the circles is proportional 
to the Li abundances quoted in Paper I. Crosses refer to 
stars without any available Li abundance. Evolutionary 
tracks, as well as the dashed and dotted lines are as in 
Fig.l 



tween 1 and 4 M Q and for different metallicities consistent 
with the range of our sample stars (see Paper I) . However, 
solar composition being relevant to most of the objects 
of our sample (about 65%), only tracks computed with 
[Fe/H]=0 will be displayed in the figures. The evolution 
was followed from the Hayashi fully convective configura- 
tion. We used the radiative opacities by Iglesias & Rogers 
(1996), completed with the atomic and molecular opaci- 
ties by Alexander & Ferguson (1994). The nuclear reac- 
tions are from Caughlan & Fowler (1988) and the screen- 
ing factors are included according to the prescription by 
Graboske et al. (1993). No transport processes except for 
the classical convective mixing (with a value of 1.6 for the 
mixing length parameter) are taken into account. 

2.3. Discrimination between dwarfs, subgiants and giants 
among the sample 

In Fig. |l| we compare the observational HR diagram 
with the evolutionary tracks computed with [Fe/H] = 0. 
The dashed line indicates the evolutionary point where 
the subgiant branch starts and which corresponds to the 
hydrogen exhaustion in the stellar central regions (i.e., 
turnoff point). About 30 stars are located below the 
turnoff line and therefore appear to be genuine dwarfs, 
although classified as subgiants in the Bright Star Cata- 
logue. On the other hand, about 15 stars located on the 
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Table 1 . Derived masses and luminosities for our program 
stars 



3.8 3.7 
Log CO 



3.6 



Fig. 3. Lithium abundances as a function of log(T e ff) for 
all the stars in our sample (Paper I). Open and filled sym- 
bols represent single and binary stars respectively. The 
circle sizes are proportional to the metal content 



right side of the dotted line have started the ascent of the 
RGB and are thus considered as giants. 



3. The Lithium main features 

Fig. |^ shows the distribution of Au on the HR diagram 
and Fig. || presents Au versus log T c ff for three ranges of 
metallicity. The main features presented in Paper I, i.e., 
the lithium discontinuity around log(T e g;) equal to 3.75, 
and the dispersion in lithium abundances for subgiants 
hotter than this value, clearly appear in both figures. 

3.1. Lithium discontinuity 

The observational lithium discontinuity actually sim- 
ply reflects the well-known dilution that occurs when the 
convective envelope starts to deepen after the turnoff and 
reaches the inner free-lithium layers (Iben 1967a, b). In 
Fig. ^ we show this behaviour in our models as a function 
of the effective temperature for different masses. A closer 
look at the models shows that the beginning of the theo- 
retical lithium dilution is a function of stellar mass, as can 
be seen in Tab. |^ where Tj^ 1 1 is the effective tempera- 
ture at the beginning of the lithium dilution for each mass. 
We also give T^ Qdu , which is the effective temperature at 
which a decrease of the surface lithium abundance by a 
factor 10 compared to the value on the main sequence is 
achieved. The dilution is a very fast process, both in terms 
of age and effective temperature interval. The theoretical 
beginning of the lithium dilution is in good agreement 



i \ i j 


l^nY TIT \ 


1V1 / 1V10 


400 


0.45±0.02 


1.0 


645 


1.36±0.05 


1.8 


3229 


0.99±0.03 


1.6 


3303 


1.45±0.04 


1.5 


4813 


0.21±0.00 


1.2 


5268 


1.71±0.08 


1.5 


6269 


1.65±0.08 


2.4 


6301 


0.64±0.03 


1.3 


9562 


0.60±0.01 


1.5 


10142 


1.69±0.06 


2.4 


11443 


1.11±0.02 


1.7 


12235 


0.56±0.03 


1.2 


12583 


1.70±0.07 


2.3 


13421 


0.92±0.03 


1.4 


16417 


0.45±0.02 


1.1 


18907 


0.65±0.02 


1.3 


19826 


1.34±0.05 


2.2 


26913 


-0.16±0.02 


1.1 


26923 


0.05±0.02 


1.2 


27536 


1.17±0.04 


2.2 


29613 


1.49±0.04 


2.2 


32503 


2.23±0.09 


3.0 


33021 


0.41±0.02 


1.1 


34411 


0.28±0.01 


1.1 


34642 


1.19±0.02 


1.8 


41700 


0.22±0.01 


1.0 


48737 


1.01±0.01 


1.5 


60532 


0.97±0.02 


1.3 


61421 


0.81±0.02 


1.3 


66011 


0.99±0.05 


1.5 


72954 


1.05±0.03 


1.6 


73752 


0.55±0.01 


1.2 


75487 


0.90±0.04 


1.4 


78154 


0.60±0.02 


1.2 


78418 


0.58±0.03 


1.2 


80956 


2.09±0.14 


3.5 


82074 


0.98±0.03 


1.5 


82210 


1.21±0.01 


1.8 


82328 


0.88±0.01 


1.4 


82543 


1.72±0.09 


2.4 


82734 


2.05±0.06 


3.5 


89010 


0.55±0.02 


1.0 


89449 


0.62±0.01 


1.3 


92588 


0.59±0.03 


1.2 


94386 


1.34±0.06 


1.3 


99491 


-0.08±0.02 


0.8 


100219 


0.61±0.03 


1.3 


102713 


1.27±0.04 


1.8 


104055 


2.18±0.14 


3.0 


104304 


-0.00i0.00 


1.0 


104307 


2.05±0.11 


2.3 


105678 


1.08±0.04 


1.6 


107295 


1.96±0.11 


3.5 


107700 


1.86±0.05 


2.6 


119992 


0.36±0.02 


1.2 


120136 


0.46±0.01 


1.2 



4 

(continued) 



J. D. do Nascimento Jr el al. : Lithium and rotation on the subgiant branch. II 
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Fig. 4. The deepening (in mass) of the convective envelope 
is shown as a function of the decreasing effective tempera- 
ture (first dredge-up) for 1.0 (solid), 1.2 (dot), 1.5 (short 
dash), 2.0 (long dash), 2.5 (dot - short dash), 3.0 (dot 
- long dash) and 4.0 M Q (short dash - long dash ) and 
[Fe/H] = 0. We present also a zoom of the region 3.80 > 
T e s > 3.74. The point labeled a indicates the end of the 
first dredge-up 

with the observed abundance drop-off along the subgiant 
branch, as can be seen in Fig. |^. 

For the considered stellar masses, the predicted dilu- 
tion factor at the end of the dredge-up (see point a on 
Fig. ^) ranges between 20 and 60 (these values obtained 
from our models do not significantly differ from the pre- 
dictions by Iben 1967a, b). This corresponds to the upper 
envelope of the observations on the right side of the dis- 
continuity, assuming a cosmic lithium abundance of 
3.1. 

3.2. Lithium dispersion 

In the following, we discuss the lithium dispersion on 

both sides of the discontinuity as a function of the stellar 
mass inferred from the evolutionary tracks at the corre- 
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Fig. 5. Lithium abundances as a function of log(T e ff ) for all our sample stars. Open and filled symbols represent single 
and binary stars respectively. The circles correspond to lithium detection while inverted triangles are for upper limits 
in the lithium abundance determination. Squares point out the main sequence stars 



Table 2. Relevant characteristics of the models at 
[Fe/H]=0 shown in Fig. ||. Column 1 gives the stellar 
masses. Columns 2 and 3 give the effective temperature at 
the point where the convective envelope just starts deep- 
ening and at the onset of the theoretical lithium dilution. 
The effective temperature at the point where the lithium 
abundance has decreased by a factor 10 compared to its 
value at the end of the main sequence is indicated in col- 
umn 4. Column 5 gives the mass of the convective envelope 
at the effective temperature of the observed rotational dis- 
continuity 
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sponding metallicity. We select four ranges of mass ac- 
cording to the ones defined in Balachandran (1995) for 
the cluster and field stars with respect to the so-called 



lithium dip region (Boesgaard & Tripicco 1986). Each sub- 
sample corresponds to peculiar observational behaviour of 
the A^i on the main sequence (see Fig. |^ where both main 
sequence and subgiants stars are shown) . In the four mass 
ranges a large dispersion of the lithium abundance appears 
and is independent of the single or binary status. 

(i) Stars with masses < 1.2 M Q show different de- 
grees of lithium depletion, which occurs already on the 
pre-main sequence and on main sequence, as observed in 
open clusters and field stars (Soderblom et al. 1993, Jones 
et al. 1999 and references therein). This explains their low 
lithium content when they reach the subgiant branch (see 
Fig. la). 

(ii) Stars with masses between 1.2 and 1.5 M Q corre- 
spond to the so-called dip region (Boesgaard & Tripicco 
1986). 

We separate the stars originating from the hot side of 
the dip region in two mass ranges : 

(iii) stars with masses between 1.5 and 2.25 M©, and 

(iv) above 2.25 M©. 

We now discuss in more detail the last three cases. 
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3.2.1. Li-dip stars 



Late F- field and open cluster dwarfs with T G ff around 
6700 K are highly lithium depleted. We find this feature 
among the dwarfs of our sample which have masses be- 
tween ~ 1.2 and 1.5 M (Fig. |b), in agreement with the 
observations in the open galactic clusters older than ~ 200 
Myr (Balachandran 1995). This lithium depletion persists 
in our subgiants of the same mass interval, in agreement 
with the observations in slightly evolved stars of M67 (Pi- 
lachowski et al. 1988, Balachandran 1995, Deliyannis et 
al. 1997). Explanations relying on the nuclear destruction 
of lithium are thus favoured by these data (see Talon & 
Charbonnel 1998 and references therein). 

3.2.2. Li in stars originating from the hot side of the dip 

As can be seen in Fig. ||c,d (see also Fig. |J), a large lithium 
dispersion exists among our most massive stars, which 
show lithium depletion by up to two orders of magnitude 
before the start of the dilution at Log T e g ~ 3.75 (point 
T e g 9 1 ). Very few observational data are available in the 
literature for stars with masses higher than 1.5 M©. In the 
Hyades, while on the main sequence these objects show 
lithium abundances close to the galactic value, except 
for a few deficient Am-stars (Boesgaard 1987, Burkhart 
& Coupry 1989). However, our observational result is in 
agreement with the findings by Balachandran (1990) and 
Burkhart & Coupry (1991) of a few slightly evolved field 
stars originating from the hot side of the dip and showing 
significant lithium depiction. We thus confirm the sug- 
gestion by Vauclair (1991, see also Charbonnel & Vau- 
clair 1992) that some extra-lithium depletion occurs inside 
these stars when they are on the main sequence, even if 
its signature does not appear at the stellar surface at the 
age of the Hyades. Among some effects suggested one can 
quote an unusual mass loss on the main sequence (Boes- 
gaard et al. 1977), and rotational induced mixing (Char- 
bonnel & Vauclair 1992, Charbonnel & Talon 1999). 

The behaviour observed in our most massive stars 
which have not yet reached the beginning of dilution ex- 
plains the very low lithium content of most of the mas- 
sive subgiants which cannot be accounted for by dilution 
alone. The underlying destruction mechanism should also 
be responsible for the lithium behaviour observed in the 
Hyades giants (Boesgaard et al. 1977) which have masses 
higher than 2.25 M Q and thus correspond to the data we 
show in Fig. |^d. This process should however preserve the 
boron abundance, which is in agreement with the stan- 
dard dilution predictions in the underabundant Li giants 
of the Hyades, as shown recently by Duncan et al. (1998). 
Observations of boron in our sample stars have to be done 
to confirm the similarity between field and cluster evolved 
stars. 

Finally, we note that Non-LTE effects alone can not 
explain the very low Li abundance derived for several mas- 
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Fig. 6. Rotation for the complete sample. The symbol 
size is proportional to the rotational velocity measure- 
ments (Vsini, in km.s" 1 ) obtained with the CORAVEL 
spectrometer by De Medeiros & Mayor (1999). The rota- 
tional discontinuity on the subgiant branch (see Paper I) 
is indicated by the two arrows. Single and binary stars are 
identified by open and filled circles respectively 

sive stars, as proposed by Duncan et al. (1998). Indeed, 
Carlsson et al. (1994) have showed that Non-LTE effects 
can reach up to ~ +0.3 dex only for such cool stars. 
This correction factor is too low to reconcile the derived Li 
abundances with the standard dilution predictions. There- 
fore, this confirms that these stars are indeed Li-poor as 
well as the hotter ones for which only upper limits have 
been derived (Non-LTE effects are much smaller for such 
hotter stars). 

4. The Vsini main features 

In Fig. 6 we present the behaviour of rotation in the HR 
diagram for our sample stars, including the dwarfs and 
giants. We define several Vsini intervals as in Lebre et al. 
(1999): Slow rotators correspond to Vsini < 10 km.s -1 , 
moderate rotators to 10 km.s -1 < Vsini < 40 km.s -1 , and 
high rotators to Vsini > 40 km.s -1 . Fig. 6 shows clearly 
the now well established rotational discontinuity along 
the subgiant branch near log(T e ff)=3.79, here indicated 
by two arrows (Gray & Nagar 1985; De Medeiros & Mayor 
1989, 1990). An interesting feature to discuss is the influ- 
ence of stellar mass on the rotational discontinuity. First 
of all, one observes that all the subgiant stars with mass 
lower than about 1.2 Mq present Vsini values lower than 
10.0 km.s -1 . Observations in young galactic clusters (see 
Gaige 1993 and references therein) show that these low- 
mass stars actually acquire a slow rotation early on the 
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main sequence, as explained by the magnetic braking sce- 
nario for main sequence stars (Kraft 1967, Schrijver and 
Pols 1993). 

Before the rotational discontinuity, the subgiants with 
mass larger than 1.2 M Q present a broad range of Vsini 
values. Because their very thin surface convective en- 
velopes are not an efficient site for magnetic field genera- 
tion via a dynamo process, these stars are not expected to 
experience significant angular momentum loss during their 
main sequence evolution. This result is again in agreement 
with the data in open clusters. 

As can be seen in Tab. || (see also Fig. [|), the effec- 
tive temperature at which the convective envelope starts 
to deepen (Tfg e ) depends slightly on the stellar mass. We 
also give the depth of the convective envelope at the ef- 
fective temperature of the rotational discontinuity, M^ ot . 
For the masses lower than 2.0 M , the observed rotation 
discontinuity occurs just when the convective envelope 
starts deepening. We thus see that if the magnetic brak- 
ing plays a relevant role in the rotational discontinuity, it 
requires only a very small change in the mass of the con- 
vective envelope. Above 2.0 M , our sample is very sparse 
(due in particular to the very rapid evolution of such stars 
in the Hertzsprung gap) and we have no data for single 
stars on the left of the rotation discontinuity exhibited by 
lower stellar masses. We cannot thus discuss further the 
impact of the deepening of the convective envelope on the 
braking in these more massive stars. 

5. Conclusions 

In this work we have analyzed the Li and rotation observa- 
tions for 120 subgiant stars from Lebre et al. (1999). We 
used the HIPPARCOS trigonometric parallax measure- 
ments to locate precisely our objects in the HR diagram 
and to determine the individual mass and evolutionary 
status for all stars in the sample. 

We have compared observed Li abundances with pre- 
dictions of Li dilution caused by the deepening of the 
convective envelope on the subgiant branch. Our mod- 
els show that the beginning of the theoretical lithium 
dilution is a function of stellar mass and coincides with 
the observational features. Stars with masses < 1.2 M Q 
show a large range in abundance before the turnoff, in- 
dicating lithium depletion in the previous phases. Stars 
with masses between 1.2 and 1.5 M Q (i.e., in the dip re- 
gion) show values in agreement with what is found in 
the open clusters. We note that many stars with masses 
higher than 1.5 Mq show lithium depletion up to two or- 
ders of magnitude before the start of the dilution at Log 
T e ff ~ 3.75. The process that depletes Li in these objects 
while on the main sequence should however preserve the 
boron abundance, which is in agreement with the standard 
dilution predictions in the underabundant Li giants of the 
Hyades, as shown recently by Duncan et al. (1998). Be- 
cause these specie burn at different depths than lithium, 



future observations of boron for our sample will provide 
powerful additional constraints and will confirm the simi- 
larity between field and cluster evolved stars. 

Our analysis confirms that low mass stars leave the 
main sequence with a low rotational rate, while more mas- 
sive stars are slowed only when reaching the subgiant 
branch. A very slight increase of the depth of the con- 
vective envelope seems to be sufficient for the magnetic 
braking to take place at this phase. Even if the decrease 
in the behavior of Lithium is nearly the rotational discon- 
tinuity, our interpretation of the observations shows that 
lithium and rotation discontinuities are independent. 
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